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ABSTRACT 

Hypervelocity stars (HVSs) escaping away from the Galactic halo are dynamical products of inter- 
actions of stars with the massive black hole(s) (MBH) in the Galactic Center (GC). They are mainly 
B-type stars with their progenitors unknown. OB stars are also populated in the GC, with many 
being hosted in a clockwise-rotating young stellar (CWS) disk within half a parsec from the MBH and 
their formation remaining puzzles. In this paper, we demonstrate that HVSs can well memorize the 
injecting directions of their progenitors using both analytical arguments and numerical simulations, 
i.e., the ejecting direction of an HVS is almost anti-parallel to the injecting direction of its progenitor. 
Therefore, the spatial distribution of HVSs maps the spatial distribution of the parent population of 
their progenitors directly. We also find that almost all the discovered HVSs are spatially consistent 
with being located on two thin disk planes. The orientation of one plane is consistent with that of the 
(inner) CWS disk, which suggests that most of the HVSs originate from the CWS disk or a previously 
existed disk- like stellar structure with an orientation similar to it. The rest of HVSs may be correlated 
with the plane of the northern arm of the mini-spiral in the GC or the plane defined by the outer 
warped part of the CWS disk. Our results not only support the GC origin of HVSs but also imply 
that the central disk (or the disk structure with a similar orientation) should persist or be frequently 
rejuvenated over the past 200 Myr, which adds a new challenge to the stellar disk formation and 
provides insights to the longstanding problem of gas fueling into massive black holes. 
Subject headings: Black hole physics — Galaxy: center — Galaxy: halo — Galaxy: kinematics and 
dynamics — Galaxy: structure 



1. INTRODUCTION 

In the past several years, surveys of hyperveloc- 
ity stars (HVSs) have found 16 HVSs with velocities 
substantially highe r than the Galacti c escape velocity 
(iBrown et all 120051: 'Hirsch et all 120051 : lEdelmann et all 
[2)05; Bro wn et al.n 2007. 20093). Most HVSs are B-type 
stars probably with masses ~ 3 — 4M0, ages > (1 — 2) x 

10^ yr and l ifetime ~ (1.6 - 3.5) x 10® yr (jBrown et all 
120071 l2009aD . Their hehocentric distances range from 
several ten to a hundred kpc and their s patial distri- 
bution o n the sky is probably anisotropic ([Brown et al.l 
I2009al[bl : lAbadi et a l." 2009). The young nature and the 
spatial anisotropic distribution of these stars, as well as 
their hypervelocities, should be related to their origin. 

A few dynamical mechanisms, involving interactions 
with the massive black hole (MBH) in the Galac- 
tic Center (GC), were propose d to eject stars with 



parameters (lYu fcTremaind 120031: IBromlev et all 120061 : 
iSesana et al.l 1200^ Hence, the spatial distribution of 



HVSs was proposed to be useful in identifying the ejec- 
tion mechanism responsible for the observed HVSs.^ 

The progenitors of the discovered B-type young HVSs 
and their sources, however, should be distinct popu- 
lations from the majority of the GC stars (which are 
typically old). Observations have indeed shown various 
young stellar structures in the GC, which may be possi- 
ble parent populations of the HVS progenitors, including 
a clockwise young stellar disk (CWS) and the other pos- 
sible c ount erclockwi se one within 0.5 pc from the central 



MBH (iLevi n & Bel oborodoA 
iLu et al.ll2009c iBartko et all 



2003t 



Paumard et"ani2006l : 



2009a), young stellar clus- 



such hyper vel ocities (Hills 19881; [Yu & Tremaind 



2006 



Bromlev et al.l 2006: Gualandris et al. 2005; Levin_ 

Baumgardt et al.. ,2006: Sesana et al.. ,2006; ,Perets et al 



2003 



20071: lo'Learv fc Loebl 120081: ILockman fc Baumgardt 



20081 ). including tidal breakup of binary star systems by 



the central MBH and three-body interactions of single 
stars with a hypothetical binary black hole (BBH) in 
the GC. The sources of (binary) stars injected into the 
vicinity of the MBH are often (implicitly) assumed to 
be isotropically distributed in previous studies, as the 
majority of the GC stars are isotropically distributed; 
and different ejection mechanisms may result in different 
spatial distributions of HVSs, depending on whether the 
central MBH is a single or a binary and relevant binary 



ters like the Arches and the Quintuplet systems at several 
ten pc away fr om the center, and the tida l remnants of 
those clusters ([Portegies Zwart et al.l 120021 ). There also 
exist other organized structures in the GC, such as, the 
circumnuclear molecular disk, the northern arm (Narm) 
and the bar components of the minispiral at a few pc 
from the center, with which young star s may b e asso- 
ciated (Yusef-Zadeh ct al. 2000; Pauma rd et al.l [2006). 
These structures may also be related to the sources of 
HVSs. All the structures above are not isotropically dis- 
tributed, but either planar or orbiting on some specific 
planes around the central MBH. The HVSs, if ejected 
by interactions of stars originating from these sources 

^ Some other observational properties of HVSs, such as the 
binarity, rotational velocity and metallicity, were also proposed 
to be usef ul in ide ntifyin g the e jection mechanisms of the 
HVSs (e.g.. ILu et al.l | 2007; Hansen 'WH; Prz vbilla eFaTI |2008| : 
ILopez-Morales fc Bonano s 2008 ; Percts 2 009a,b,) . 
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with the central MBH(s), are very likely to be spatially 
correlated. And their spatial distribution should be map- 
ping the distribution of the parent populations of their 
progenitors, if the ejected HVSs can well memorize the 
injecting direction of their progenitors. 

This paper is organized as following. In § [21 we study 
how the spatial distribution of HVSs is related with the 
geometrical structure of the parent population of their 
progenitors. In § [3] we find that most HVSs discovered 
so far are probably originated from the clockwise young 
stellar disk (CWS) or a disk-like stellar structure with an 
orientation similar to that of the CWS disk in the GC. 
Discussion and conclusions are given in § 3] and § O 

2. THE DEFLECTION ANGLE: THE DIRECTION CHANGE 
OF HVSS FROM ITS PROGENITOR 

In this section, we use both analytical analysis and nu- 
merical calculations to demonstrate that HVSs can well 
memorize the injecting direction of their progenitors. We 
do this for different possible dynamical mechanisms of 
ejecting HVSs introduced in § (T] 

2.1. Tidal breakup of a binary star by the central MBH 

If a binary star initially unbound or weakly bound to 
the central MBH approaches the vicinity of the MBH 



within a tidal distance Rti 



•tid 



ah 



( 3M. \ 

\niHVS+™c J 



1/3 



the bi- 



nary is probably tidally broken up, and one component 
of the binary may be ejected as a HVS, where M, is 
the m as s of the central MBH [~ 4 x IO^Mq, iGhez et al.l 
(|2008f) : iGillessen et"all (|2009[ )]. torvs and rUc are the 
mass of the ejected HVS and the other binary compo- 
nent, r espectively, and Oh is the semi-ra ajor axis of the 
binary (|Hillslll988i : [Yu fc Tremaindl2003f l. The initial in- 
jecting velocity of the binary is <C \/GM,/Rtid, where G 
is the gravitational constant, and then t he velocity of the 
ejected HVS can be appr oximated by (|Yu fc Tremaind 
12003; .Bromlev et al.l[200l 



?;?ys ^ km s" 
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(1) 



The orbit of the ejected HVS and the initial orbit of the 
injecting binary both have eccentricities e close to 1 and 
can be approximated as rectilinear at distances faraway 
from the MBH. Below we aim to find out the range of 
the deflection angle of these two rectilinear directions. 

For convenience, we first consider a purely two-body 
problem that a star with mass r7i»(^ Af,) starts at a 
velocity — {GM,/aY^'^ from infinity and is mov- 
ing towards the MBH on a hyperbolic orbit with pe- 
riapsis distance i?min o,. The velocity of the star 
at periapsis is Vp = Voo{2a/ Rmin + 1)^^^, and then the 
deflection angle of its direction from its initial inject- 
ing velocity is tt/2 — ^Pi, where tan^*! = ^^jl — e^| ~ 
y/2R~J^ ~ w™(2i?„iin/GAf,)i/2. If the star moves to- 
wards the MBH on a parabolic or elliptical orbit with 
eccentricity close to 1, the deflection angle of the stellar 
velocity moving from almost the infinity (excluding loca- 
tions near the apoapsis) to the periapsis i?min is about 



tt/2. Similarly as above, for a star that can escape the 
BH with velocity v^^, the deflection angle of the es- 
caping star from its velocity at i?niin is tt/2 — ^2 with 
tank's ~ v^^{2R^in/GM,y/'^. For the HVSs discov- 
ered so far, we have ~ 750 — 1000 km s~^ obtained 
by removing the velocity deceleration due t o the Galac- 
tic potential measured bv iXue et al.l ()2008D . Note that 
the relative change of the eccentricity vector (pointing 
towards the periapsis from the MBH) is ^ Svp/vp ^ 
and it is negligible compared to ^'2 for sufficiently small 
*2, where Svp - { mHvs+rnJ VG{m^^s +mc)/ab is the 
change of the velocity of the binary component ejected 
as the HVS after the binary breakup. For the process 
of tidal breakup of a binary star at a distance of i?tid 
from the MBH, we approximately have i?min ~ -Rtidj and 
the total deflection angle 8 of the ejected HVS from the 
original injecting binary is about tt — y^'^i + ^i- 
9 is ~ TT - (^'i -I- \E'2) if the HVS and the injecting bi- 
nary are on the same orbital plane and tt — ^2 for 
ug^vs ^ ^;ini^ Pqj. _ 0.6 AU, wg^s ^ 1000 km s^i, and 

ug^* ~ 250 km s~\ we have ^Z*? + *i 0.2 - 10°, that 
is, the HVSs are almost reversing the injecting direction 
of their progenitors. 

To confirm the above analysis, we numerically real- 
ize the process of ejecting an HVS as a binary inter- 
acts with a MBH. We use an explicit 5(4)-order Runge- 
Kutt a scheme to integrate the full three-body prob- 
lem (jPormand fc Princd llQSGl : iHaier et all I1993D . In 
the three-body simulations, we first assume that all 
binary stars initially move on hyperbolic orbits with 
w™' ~ 250 km from infinity. In the calculations we 
set other relevant parameters as follows: (1) The dis- 
tribution of the semi-major axes ab of binary stars is 
assumed to follow P{ab)dab (x -^dab as suggested by 
obser vations of binaries with 0- type or B-type primary 
stars (|Kobulnickv fc Frverl l2007n . The lower limit of ab 
is roughly set to 0.03 AU, i.e., about twice the physi- 
cal radius of a 4AfQ star; and the upper limit of Of, is 
set to 2 AU to ensure the ejected star can escape to 
large Galactic radii. (2) The mass distribution of pri- 
mary stars mi follows the Miller-Scalo in itial m ass func- 
tion, i.e., /(mi) cx m||~" and a ~ 2.7 (jKrou pa 2002). 
For massive binary stars, the distribution of the sec- 
ondary star or the mass ratio q = m2/mi(< 1) can 
be described by two populations: (a) a twin population, 
i.e., about 40% binary stars have g ~ 1 or m2 ~ mi, 
and (b) the rest binaries foll ow a /(q) ^ constant distri- 
bution C Kobulnickv fc Frveri [20071 : IKiminki et al.l 120081 . 
2009). (3) The eccentricity of the binary star is as- 
sumed to be 0. (4) The orientation of the binary orbital 
plane is randomly chosen. (5) The probability distri- 
bution of the closest approach distances to the MBH is 
^'(-Rmin)'^-Rmin oc c?i?niin, which corrcspouds to the im- 
pact parameter distribution p(b)db oc bdb. We only se- 
lect those cases that the masses of ejected stars are in 
the mass range (3AfQ,4Af0) of the observed HVSs and 
calculate the probability distribution of their deflection 
angles 9 . As seen from panel a of Figure [1] the values of 
9 range from 165° to 180°, consistent with our analysis 
above. The distribution of 9 does not show dependence 
on . The reason for this independence is as follows. 
If the mass of the central MBH (here the MBH in the 
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GC) and the mass of an HVS are fixed, ug^^ is primar- 
ily determined by the semi-major axis of the binary, and 
given i?tid oc Ob, tan ^2 oc 

(2i?tid/G'M.)i/2 is inde- 
pendent of Ob and hence . 

We also calculate the deflection angle distribution for 
the case that the injected binary stars are initially on 
weakly bound and highly eccentric orbits, instead of un- 
bound orbits. In this case, the only difference in the ini- 
tial conditions from that for initially unbound binaries 
is as follows: the apoapsis distribution of the orbit of 
the binary barycenter follows p(rapo)d?'apo oc r'^^^dr^^po, 
Tapo is in the range (0.04 pc, 0.5 pc), and the distri- 
bution of the closest approach -Rmin, i-e., the periapsis 
distance of the orbit of the binary barycenter, is the 
same as above. The above distribution of the apoapsis 
is adopted so that the energy distribution of the binaries 
injected from a disk structure is consistent with the den- 
sity distribu tion of young stellar disk rec ently discovered 
m th e GC IL evin & Bcloborodov 2003: jPaumard et al.l 
[2OO6L iLu et al 1 120091 : iBartko et all i2009"ah . The radial 
distribution of sta rs within t he yo ung stellar disk plane 
is g r-'^-3±o.7 i n |Lu et all (gOOl) and oc r'^-^^^-^ in 
iPaumard et al.l ()2006[ ). The resulted distribution of 9 
is shown in panel b of Figure [TJ And 8 is also in the 
range from 165° to 180°, which again confirms our anal- 
ysis above that HVSs memorize the direction of their 
original binaries well. 

2.2. Dynamical ejection by a binary black hole 

HVSs can also be produced by three-body interactions 
of a single star with a hard BBH. A BBH is hard if 
its semi-major axis obbh is less than ah — G'M,^2/4cr^, 
where a is the stellar velocity dispersion of the host 
galaxy and M,^2 is the mass of the secondary black hole. 
For a hard BBH, most low-angular-momentum stars that 
can enter into the region r < obbh will be ejected after 
one or several encounters with the BBH and the r.m.s. 
of the velocities of the ejected stars at infinity is 



2KGM,,iM.,2/iM.aBBii) 



930 km s-'m°-25(l - zy)^/"(0.1ah/aBBH)'^',(2) 



1/2/ 



.1/2 



[see eq. 1 in iLu et all ()2007f) ]. where Af.^i is the mass of 
the primary black hole, M, = A/,^i+A/,^2, = M,,2/M,, 
m, = M,/(4 X IO^Mq), ~ 1.6. Unless the eccentric- 
ity of the BBH is excited to an extr emely high value 
shown in some numerical simulations (|Baumgardt et al" " 



20061: iMatsubavashi et all 120071 : iLockman fc Baumgardt 



20081) , the BBH can s tay at its hard stage fo r a long time 
(e.g., up to 10^ yr in lYu fc Tremain"il2003[ ) and ejection 
of HVSs from the GC would last that long. A BBH with 
an extremely high eccentricity has a much shorter life- 
time due to gravitation radiation, which is substantially 
smaller than the travel time span (^^ 2 x 10* yr) of the 
observed HVSs. 

In this mechanism, HVSs are ejected out from dis- 
tances i?min ^ obbh, and we have tan\I'2 — 



2fi„ 



GM. 



< 2.5 



If M,_2 ~ we have 

\l/2 ~ 50°, and thus HVSs have lost their memory of 
the direction of their progenitors and the information of 
the BBH orbital plane may be imprinted on the spa- 
tial distribution of HVSs. Current observational con- 



straints on the mass of the secondary BH in the GC 
gives M., 2 < O.OIM.,1 (jHansen fc Milosavlie"^ 120031 : 
lYu fc Tremainell2003t lGillessen et al."2009^. for which we 
have ^'2 < 14°. If the progenitor of the HVS is unbound 
to the BBH with small initial velocities, it is ejected out 
generally after one or a few close encounters with the 
BBH, and we have tan^-i < w™(aBBH/GA/.)i/2. For 

<C , we have <C VE'2- If the injecting star 
is initially on a bound (and highly eccentric) orbit, it 
may be ejected out as a HVS after many times of en- 
counters with the BBH, and the accumulated relative 
change of the eccentricity vector is ^ ^'2 and negligible 
for sufficiently small '52. In this case the total deflection 
angle 8 of the ejected HVS from the injecting direction 
of its progenitor is thus tt — ^J^>\ + ^>\ ^ it — '^2, and 
thus HVSs may have a good memory of the direction 
of their progenitors. We do numerical experiments on 
three-body interactions between single stars and a BBH 
with M. = 4 X IO^Mq, v = 0.003, and obbh = O.lah, 
where the two BHs are set to be sufficiently close so that 
they are able to eject stars with the observed hyperve- 
locities. As shown in Figure [21 the defiection angles 8 
obtained from the calculations range from 165° to 180°, 
no matter whether the injecting stars are initially un- 
bound (panel a in Fig. [2]) or weakly bound (i.e., highly 
eccentric) orbits (panel b in Fig. [2|). The result is not 
sensitive to the BBH eccentricity and abbh- 

Hypervelocity stars may also be produced by interac- 
tions of a single star with a stellar-mass black hole in 
the vicinity of the central MBH. In this case the veloc- 
ity kick 5v received by the star when it passes by the 
stellar-mass black hole is typically larger and compara- 
ble to the orbital velocity of the stellar mass black hole 
Uorb and thus the deflection angle deviates from tt sub- 
stantially larger than that for the cases discussed above. 
Therefore, it may be difficult for the HVSs produced by 
this mechanism to memorize the direction of their pro- 
genitors and hence difficult to explain the consistence of 
the HVS plane and the CWS plane shown in § [3] below. 

Additional deflection of the moving direction of HVSs 
may be introduced if the Galactic potential is substan- 
tially flattened or triaxial; but this deflect ion, typically 
< 2° for v^'' ^ 1000 km s-\ is negligible (|Yu fc Madaul 
IMiTHG^in et al.l[20f)5h . 

In summary, it is plausible to use the spatial distri- 
bution of HVSs unbound to the Galactic potential to 
map the parent population of their progenitors and re- 
veal their origin if the HVSs are produced by the tidal 
breakup of binary stars or the BBH mechanism. 

3. COMPARISON WITH OBSERVATIONS 

Figure [3] shows the spatial distribution of all (16) HVSs 
detected so far {filled circles) in the Galactic coordinates 
by a Hammer- Aitoff projection. The spatia l distribution 
of thes e HVSs is anisotropic at a 3.5-(t level dB rown et aLl 
|2009ai rbl) and the majority of the H VSs are located a t 
Galactic longitudes I - 240° - 270° (jAbadi et al.ll2009( ). 
The projection of other various planar structures in the 
GC are also plotted in the figure. As shown in Figure [3l 
most HVSs (11 among 16) situate close to the plane of 
the clockwise young stellar disk; and a second group, 
including 4 other HVSs and possibly a specific one among 
the above eleven, situate close to the Narm plane. Only 
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Fig. 1. — Probability distribution of the deflection angle © of 
HVSs for the mechanism of tidal breakup of binary stars. Panel a 
is for the case that the injected binary stars are initially unbound 
to the MBH, panel b for the case that the i njected binary stars are 
initially on weakly bound orbits (see 12.11 1. The solid and dashed 
lines represent the distribution of deflection angle for HVSs with 
velocities in the range of (750 km s~^,1000 km s~^) and 

(1000 km s~^, 1500 km s~-^), respectively. 



0.25 




Fig. 2. — Probability distribution of the deflection angle © for 
the mechanism of three-body interactions between single stars and 
a hypothetical BBH in the GC. The semi-major axis of the BBH 
is O.loh and its mass ratio is 0.003. Panel a is for the case that the 
injected single stars are initially unbound to the central BBH, and 
panel b for the case that the injected stars are initially on weakly 
bound orbits. Different line types have the same meaning as those 
in Fig. 1. Thick lines are for a circular BBH and thin lines for a 
highly eccentric (e = 0.9) BBH. 

one HVS, i.e., HE 0437-5439, is neither on the CWS plane 
nor on the Narr n plane, which was sum ested to be ejected 
from the LMC (jEdelmann et al.l[2005f ). Almost all HVSs 
are quite far away from the planes of other structures 
plotted in the figure. 

We use the statistic to fit the possible planar struc- 
tures of the HVSs. The normal n to the best-fit plane of 

the HVSs is obtained by minimizing = ^i( ft-ds^ )^ 
(i ~ 1,2...), where eV- is the unit position vector of 
the i-th HVS seen from the GC and dcn is its possi- 
ble deviation from an original plane for which we take 
equivalently as the error of eV. in standard statistics. As 
the thickne ss (half-opening angle) of the CWS plane is 
~ 7° - 10° (|Lu et al.ll200a iBartko et al.ll2009aD and the 
deviation angles of ejected HVSs from the direction of 
their parent population should be about ^ 0° — 15°, we 
set the amount of de^i correspondingly to a deviation of 
7°. Thus, we obtain the normal of the best-fit plane of 
the eleven HVSs (/,6) = (309°, -15°) with = 0.82, 
where P^2 gives the probability of being higher by 
chance and it is high enough so that the fit is accept- 
able. We take the error of the fitted normal direction as 
61 ±6° and 6b ~ ±8°, which roughly correspond to 
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Fig. 3. — Hammer-Aitoff projection in the Galactic coordinates 
for all the discovered HVSs and some organized structures in the 
GC. Open red circles represent the coordinates (centered on the 
sun) of these observed HVSs, while flUed red circles represent their 
coordinates if projecting the HVS positions to inflnity on the sky 
of an observer sitting in the GC. The curves show the planes (also 
projected to inflnity) of the clockwise young stellar disk (CWS; 
magenta), the counterclockwise disk (yellow), the northern arm 
(Narm; blue) and the bar (cyan) components of the minispiral, and 
the circumnuclear disk (green), respectively (Paumard et al. 2006). 
The magenta curves represent the fitted plane of the clockwise 
young st ellar disk from Paumard et al. (2006) (thick solid curve), 
ILu et aU t2009. ) (thin dashed curve) and .Bartko et al., (i2009a ) (thin 
dotted curve), respectively. 
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Fig. 4. — Similar to Fig. |3] with adding the four HVS candi- 
dates (carmine stars), the eight bound HVS (cyan stars), and the 
plane of the outer warped part of the CWS disk recently reported 
by Bartko et al. (2009a) (yellow dashed curve). For clarity, some 
curves in Fig. [3] are removed in this figure. See details in §|5] 

the 68% confidence level. The normal of the CWS disk 
plane, or at least the normal of the CWS disk in the inner 
regio n (0.8"-3.5"), (/ 6) = (310°. -18°V. jPaumard et all 
2006; ILu et "aI1l2009l : iBartko et al.ll2009aD are within the 
error range. For the secondary group of five HVSs, 
we obtain (l,h) = (188°, -52°) with P^^ = 0.87, and 
{61, 6b) = (±22°, ±5°) correspond to the 68% confi- 
dence level. The no rmal to the observed Narm plane 
(/, b) = (162°, -47°) (jPaumard et al. 2006) is at the 80% 
c onfidence leve l of the best-fit value. 

IBartko et al.l ()2009a[ ) recently reported that the outer 
part of the CWS disk may be significantly warped and 
the normal of the warped part is {l,b) = (136°,— 44°). 
This normal is close to the normal of the Narm plane as 
shown in Figure [H which suggests some physical connec- 
tion between the Narm and the outer warped part of the 
CWS disk. The normal of the outer warped part of the 
CWS disk is at the 99.9% confidence level of the best fit. 
Compared with the plane of the Narm, statistically the 
plane defined by the normal of the warped part of the 
CWS disk is less likely to be consistent with the fitted 
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plane of the second population of HVSs. 

It is worthy to n ote that there are f our other HVS 
candidates listed in 'Brown et al.l (|2009af ). As shown in 
Figure H two of them (SDSS J1403+1450 and SDSS 
J1546+2437) are within - 16° to the CWS plane and the 
other two (SDSS J0940+5309 and SDSS J1014+5631) 
are within 3° to the Narm plane. If they are confirmed 
to be HVSs and included in the fit, the normals to the two 
best-fit planes are {l,b) = (318°, -9°) with P^2 = 0.81 
and {l,b) = (180°, -50°) with P^2 = 0.79, consistent 
with the planar structures fitted above. Even for the 
eight b ound "hypervelocity" stars listed in iBrown et al.l 
(|2009af) . four of them are close to the CWS disk within 
15°, so are three of them to the Narm plane (see Fig. 2]). 
After including them, the normals to the two best-fit 
planes are {l,b) = (311°, -14°) with P^2 = 0.44 and 
{l,b) = (176°, -53°) with P^2 = 0.30, respectively. As 
seen from Figure [51 after including the HVS candidates 
and the bound ones, the planar structure close to the 
CWS disk appears more obvious as their locations ex- 
tend to a wider area in the sky. Only one object among 
the bound sample, SDSS J1404-I-3522, with the smallest 
Galactocentric distance, is significantly separate from the 
above two planes. We note here that the bound popula- 
tion of "hypervelocity" stars are more likely to be con- 
taminated by the high-velocity stars produced by other 
mechanisms as they have substantially smaller velocities 
than the unbound stars. 

We have also tested that the fits cannot be passed sta- 
tistically if choosing other observed structures, such as, 
the bar, the circumnuclear molecular disk or the counter 
clockwise-rotating disk. The fit cannot be passed, either, 
by fitting all the stars to one best orientation. 

Assuming that the parent population of the observed 
HVSs are on the two fitted disks with a thickness of 7° 
and normal of {l,b) = (311°, -14°) and (176°, -53°), 
we simulate the process of tidal breakup of binary stars 
around the central MBH using similar initial conditions 
as that in § [2] As shown in Figure [5l the spatial dis- 
tribution of the simulated HVSs can well match the dis- 
tribution of observed HVSs. Similar spatial distribution 
can also be reproduced if the BBH ejection mechanism is 
alternatively adopted. For simplicity, we do not present 
that in details. 

4. DISCUSSION 

As we have demonstrated above, the spatial distribu- 
tion of the discovered HVSs is consistent with being lo- 
cated on two thin disk planes, and these two planes are 
consistent with that of the (inner) CWS disk, and the 
Narm or the outer warped part of the CWS disk, respec- 
tively. We discuss two possible explanations to these 
results below. 

One explanation could be that the HVSs are originated 
from some unknown and previously existed disk-like stel- 
lar structures with orientation similar to that of the CWS 
disk, and the Narm or the outer warped part of the CWS 
disk. If this is true, one needs to answer what and where 
the unknown structures are, why their structures are con- 
sistent with the CWS and Narm (or the warped outer 
part of the CWS disk), and whether the consistency is 
coincident or some natural outcome. 

The other explanation is that most HVSs are origi- 
nated from the CWS disk and a second population of 



90° 
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Fig. 5. — The spatial distribution of simulated HVSs and the 
spatial distribution of observed HVSs. Red solid circles, carmine 
stars and the cyan stars have the same meanings as those in Fig|4] 
Open circles and diamonds represent the simulated HVSs origi- 
nated from the two fitted planes with a thickness of 7° , and normal 
{l,b) = (311°, -14°) and (176°, -53°), respectively (see details in 
§ O . The magenta and blue curves represent these two best fit- 
ted planes. This figure shows that the spatial distribution of the 
simulated HVSs matches that of the observed ones quite well. 

HVSs may be originated from the Narm or the outer 
warped part of the CWS disk. For this explanation, 
the young CWS disk should persist or be frequently re- 
juvenated over the past 2 x 10® yr as constrained 
by the travel time of these HVSs (1 - 2) x 10* yr], 
which is extremely puzzling. The ages of the OB stars 
in the disk are only ~ 6 ± 2 Myr. As the in-situ 
formation of the CWS stellar disk is already difficult 
due to the suppression of star formatio n by the MBH 
tidal field (Levin & Bcloborodov 2003; Lu et al.l 120091: 
'Bartko et al. 2009a; Sanders ,19921 : iNavakshin fc Cuadral 
2005; BonncU & Rice 2008), how have young stars been 
continuously forming in the disk and in the meantime 
how can the disk plane maintain its direction? One key 
to solve this puzzle would be continuous sinking of cold 
gas onto the CWS disk. However, the observed gaseous 
structures located just outside the CWS disk generally 
do not have the same direction as the disk. Solutions to 
this puzzle will provide profound insights to the gas fu- 
eling into the vicinity of the central MBH in the GC and 
MBHs in other galactic nuclei in general. Detailed stud- 
ies of physical properties of HVSs (e.g., metallicity) and 
their relations with those in the CWS stellar disk may 
also help to unveil the secret of their star formation. 

It appears that the counter-clockwise disk is not cor- 
related with the currently discovered HVSs. The current 
observed B-types in the CWS dis k region are a l so mor e 
isotropic than the 0/WR stars (|Bartko et al.l l2009bf) . 
which appears not to be the same as the planar distri- 
bution of HVSs. The correlation between the observed 
HVSs and the CWS disk plane may suggest that (1) only 
the B-type stars on the disk plane can be perturbed to 
inject into the immediate vicinity of the central MBH as 
the prog enitors of HVSs (e.g., by secular evolution of the 
disk, see lMadigan et al] 120091 ): and (2) B-type stars ini- 
tially formed on the disk may be heated up by relaxation 
processes later. Detailed study of this relation may pro- 
vide some constraints/hints on the mechanism to deliver 
stars to the vicinity of the central MBH. 

It is worthy to further explore how the orbits of (bi- 
nary) stars in the CWS and the young stellar structure 
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associated with the Narm plane or the outer warped 
part of the CWS disk are perturbed so that they can 
move to the immediate vicinity of the central MBH(s). 
This pe rturbation may be due to some massive per- 
turbers (jPerets et alJ l2007fl or secular evol ution of the 
structures themselves (iMadigan et all 12009( 1. Perturba- 
tions on other young stellar structures, such as the tidal 
streams of young stellar clusters like the Arches and the 
Quintuplet systems, may also inject (binary) stars to the 
immediate vicinity of the central MBH(s) and lead to 
ejection of HVSs. Therefore, there may be other planar- 
like spatial distribution of young HVSs. If these could 
be found in future HVS surveys, together with those 
found so far, the spatial distribution of HVSs will be 
mapping young stellar structures ever existed in the GC 
over the past 2 x 10* yr. With more and more HVSs to 
be discovered in the all-sky survey in the future, the sta- 
tistical methods on how to extract the stellar structures 
would need to be improved and the detailed improvement 
method should depend on how complicated or simple the 
structures would be. 

One natural prediction of the HVS origination pro- 
posed above is that some B-typc HVSs exist close to the 
CWS disk plane and the Narm plane or the outer warped 
part of the CWS disk in the southern hemisphere, which 
should be a crucial check by future HVS surveys. 

If old-population HVSs can be detected, their spatial 
distribution may be different from those shown in Fig- 
ure [3l as the parent population of their progenitors may 
be significantly isotropic than that of the discovered B- 
type HVSs. Studying the spatial distribution of different 
types of HVSs may help to reveal information on the star 
formation and the dynamical environment in the GC. 

5. CONCLUSIONS 
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